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Abstract 

In human society, a lot of social phenomena can be concluded into a math¬ 
ematical problem called the bipartite matching, one of the most well known 
model is the marriage problem proposed by Gale and Shapley. In this arti¬ 
cle, we try to hnd out some intrinsic properties of the ground state of this 
model and thus gain more insights and ideas about the matching problem. 
We apply Kuhn-Munkres Algorithm to hnd out the numerical ground state 
solution of the system. The simulation result proves the previous theoretical 
analysis using replica method. In the result, we also hnd out the amount 
of blocking pairs which can be regarded as a representative of the system 
stability. Furthermore, we discover that the connectivity in the bipartite 
matching problem has a great impact on the stability of the ground state, 
and the system will become more unstable if there were more connections 
between men and women. 

Keywords: bipartite matching, marriage problem, optimization, blocking 
pair, K-M algorithm 


1. Introduction 

Bipartite matching problem, which is to match disjoint two groups of 
agents, is widely seen in human society, for instance, the marriage prob¬ 
lem between men and women, college admission problem between students 
and universities, assignment between workers and jobs, and also the choice 
making between buyers and sellers. 
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Gale and Shapley first introduced the stable marriage problem, a one-to- 
one two side matching [l|, which is still one of the most important bipartite 
matching problems. Gale-Shapley Algorithm guarantees the existence of a 
stable solution that is also the optimal stable solution for the active side. Due 
to its multiple applications in society and its intriguing fascinating proper¬ 
ties, not only economists but also the statistic physicists are attracted by 
the bipartite matching problem. By giving up the restriction to stabilize the 
system, Parisi and his colleagues used the replica method of spin glass theory 
to study the global optimal solution of one-to-one two side matching prob¬ 
lem (3, m. Later, Zhang et al. studied the scaling behavior of Gale-Shapley 
model, partial information matching and analyzed the common features of 
all the stable solutions . Moreover, Dzierzawa introduced the accratance 
threshold and thus improved the matching result of the passive side j7|. Re¬ 
cently, Zhou et al studied the bidirectional selection problem, provided us a 
new approach from human society network 0,[i3. 

This article is organized in the following way. First we introduce some 
fundamental background about Gale-Shapley model. Then we use numerical 
simulations to study the optimal solution based on the bipartite marriage 


problem. By bringing the Kuhn-Munkres Algorithm [ll|, ll2| into the mar¬ 


riage problem, we capture the ground state of the system. Next we analyze 
the average energy and other properties of this solution in detail. After that, 
the number of blocking pairs in the ground state proposed by Zhang 0] is 
given by the simulation result, the dynamics is also discussed here. At last, 
we make a simple a quantity analysis of the stability of system ground state. 


2. Model 

Here we study the optimization solution of the bipartite matching prob¬ 
lem on the basis of the Gale and Shapley marriage model. This model con¬ 
sists of two sets of agent, ^ = {mi, m 2 ,..., m„} standing for n men and 
W = (tCi, W 2 , ■ ■ ■, iVn} for n women. A hnal outcome of the marriage prob¬ 
lem is a one-to-one matching of men and women, i.e. an invertible bijection 
X : ^ —)■ W. An outcome x can be denoted as: 


X = [(mi,x(mi)), (m 2 ,a;(m 2 )),..., (m„, a:(m„))]. 


which x{mi) = Wa means the woman who matched with man m,, and 
x~^{wa) = mi is the man matched with Wa 13 . 
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One of the most attractive questions all along these years is the optimal 
solution of the system. In order to hnd out the optimal solution, we need a 
measurement for the satisfaction of each man and woman. Compared with 
the way of nature that it always looks for the ground state, we also assign an 
energy term Si to agent z, to represent how satisfactory he/she is. The energy 
term actually represents the ranking of the assigned mate for agent i. The 
smaller e, the happier he/she is. Several previous works (ij-Q postulated the 
energy term e follows a discrete uniform distribution like e = 1,2,...,n in 
their model, but we think this might be extended to some other distributions. 
Since we have a large population 2n in the system, there would be no much 
difference between the expectation values of discrete distribution and the 
continuous distribution. So here we suggest a more soften assumption that 
the energy e is randomly uniform distributed on [0,1], and we predict that 
it will not change the result much from the previous work With the 

description of satisfaction of each agent, the optimal solution problem will 
be rephrased into a minimum energy problem. 

For the convenience, we denote the men with Latin alphabet, and women 
with Greek alphabet. We define matrix M and W, in which Mi^a is the 
energy of man rrii if women Wa was matched to him, all the same, Wpj 
is the energy of woman wp if man rrij was matched to her. For a given 
assignment x = [(mi, a;(mi)), (m 2 , x(m 2 )),..., (m^, a:(m„))], we can calculate 
the average energy of each man and woman: 
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Then we have the average energy per person: 
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Here we define the matrix H = where the element „ is the mean 

energy of the man rrii and woman Wa- With the matrix H, we now are able to 
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solve the optimal solution e™” i.e. the minimum of Eh, and the corresponding 
matching x is called the ground state. 

In the example of a matching problem which consists n men and n women, 
there will be n\ different states of the system in total. Obviously it’s neither 
possible nor clever to find out theground state by Exhaustive method. Facing 
with a similar situation, Kuhn [ll| and Munkres [l^ developed an algorithm 
called Kuhn —Munkres algorithm, to solve minimum solution in a weighted 
bipartite matching. We think this method can also be applied to our problem 
here. The total energy matrix H can be regarded as the cost matrix in 
the bipartite assignment problem. In this way, it is possible to dig out the 
matching corresponding to the ground state which has the lowest global 
energy only spending time of 0{n'^). But with previous study of the matching 
problem, it’s unfortunate that the ground state of this matching problem is 
not stationary state, which means it is not stable. Here we call a matching x 
is stable only if there was no pairs of man and woman who prefer each other 
than their assigned mate in x, that is to say, for each individual there is no 
any better choice to improve their personal situation. But if such pair of man 
ruj and woman Wa exists, we name it a blocking pair jl, 14|, abbreviated as 
BP hereafter. The more BPs a system contains, the more unstable it is. 
We are curious about the quantity of HFs in a matching problem, becausejt 
can be treated as a representative for the stability of a certain matching 


3. K-M Algorithm 

Before we give a description of this algorithm Hlli , we’d like to intro¬ 
duce two lemmas: 

Lemma 1. The equivalent statement of Konig’s theorem: in matrix H, 
the biggest amount of independent zero elements we can find is equal to the 
minimum number of lines (rows or columns) that covers all the zero elements. 

Lemma 2. Obviously, the minimal assignment solution does not change, 
if we add or subtract the same number in any rows or columns at the same 
time. 
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Algorithm description 

Step 1. Starting with a given matrix H, for each row and colnmn, 
we snbtract its minimnm valne from the rows and colnmns (thus we 
create one or several zeros in those rows and columns). We denote the 
new matrix Hi. 

Step 2. Find the smallest set of lines (vertical and horizontal), 
which consists of Uk lines, to cover all the zeros in the matrix Hi. 
Step 3. If Ufc = n, then we will have n independent zeros, the corre¬ 
sponding n positions represent the dual assignment. 

Step 4. If Uk < n, hnd out the smallest matrix element which was 
not covered by any of Lk. We add ak to all the uncovered rows and 
subtract ak from all the uncovered columns. Then return to step 2. 


4. Simulation and Analysis 


We applied this Kuhn-Munkres Algorithm to our model, and studied 
the optimal solution he.ground state of the marriage problem where the 
population of men or women n = 100,200, 300,..., 1000. 

With each fixed population of n men and n women, we did 100 trials 
starting with randomly generated H matrix, the result of the average energy 
per person in the ground state (e™"') and its standard deviation a are shown 
in Tab.l. We see that the standard deviation can be estimited as the recip¬ 
rocal of n. That is because each agent is identical to another who are weakly 
correlated, by the Central Limit Theorem 15(], we can know the e™"' follows 
the normal distribution. Furthermore, we make 10000 trials when n=100, 
the simulation result proves our analysis, shown in Fig.l. 


Table 1: Mean and Standard Deviation of Energy Per Person 


n 

too 

200 

300 

400 

500 

600 

700 

800 

900 

/-mtn\ 
\^H / 

(7 

0.0802 

0.0034 

0.0569 

0.0017 

0.0464 

0.0011 

0.0405 

0.00091 

0.0361 

0.00070 

0.0329 

0.00056 

0.0305 

0.00047 

0.0285 

0.00042 

0.0269 

0.00038 


The (e™"") for the 100 trials is shown in Fig.2. With a power exponent 
£t, we get the power of —0.495 and also the expectation of average energy 
per person is which is very close to the predicted value of Zhang 
et al with replica method. 
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Figure 1: The distribution of 10000 trials on when n = 100. 

Furthermore, the distribution of personal energy e is studied here as well, 
shown in Fig.2. By taking the semi-logarithmic, we hnd that the probability 
distribution function f{e) decays faster than exponent, but it can be htted 
almost perfectly by an exponent function with its index being quadratic 
function: 


— g-a(v/n£)^-6(^£)+c^ 

y/n 

If we treat the y/ne as the argument and the as the dependent variable 
(in this way we normalized the scaling factor), and then fit the parameter a, 
b and c in Eq. (|4]) with the simulation results, we get a = 0.19, 6 = 0.68, and 
c = 0.02, so the distribution function can be written as: 


/(g) ^ . g-0.19(Vne)2-0.68(Vne)+0.02 

\/n 

To ensure that these parameters are reasonable, we integrate the distribution 
function and it shows that: 


/(e) de = 1.03, 


( 6 ) 
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Figure 2: The relationship between (e^™) and the system scale n, for n = 
100, 200, 300,..., 1000. The dashed line indicate the fitting of 


ef{e) de = 


0.84 


n 


(7) 


The result is a little bigger than the actual value, which means the actual 

2 

distribution function decays even faster than e~‘^ , but this function is still 
in good accordance with the simulation result. 


5. Blocking Pairs 

As Zhang et al mentioned , the quantity of blocking pairs can be used as 
a measurement of the stability of a solution, and he gave a rough estimation 
of the number of blocking pairs being {neY ~ 0.65n . 

Now we try to think more about this. First let us consider the pos¬ 
sibility that a man cannot find a BP. Imagine a man i with energy e, 
there would be ne women who are possible to form a. BP with man z, 
i.e.Fi^ai-, Fi^a 2 i ■ ■ ■) Fi,ane < Fi^x{i)- If these ne women satisfied the follow¬ 
ing situation that their current mate were all better than man i, which can 
be represented as: 


there would not be any blocking pair forming with man i. 
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Figure 3: (Color online) The PDF of e with n = 100(black circles), 300(red circles), 
1000(blue circles). The dashed lines are the corresponding fitting by Eq. ([5|). 


In the process of chasing gronnd state, if we neglect other influencing 
factors and fix the Ga^i, we can see that the possibility of woman a matching 
to man i gets bigger when Fj q, gets smaller. So in the situation that woman 
a was not matched to man i even if Fj „ was very small, it’s reasonable to say 
that the expected energy of men i is bigger than average. Then the possibility 
of woman a and man i being a BP is smaller than Ga^x-^(Q)- Therefore, the 
possibility Fo(e) that man i with energy £ cannot find a blocking pair: 

ns 

-Po(e) > ]^(1 - Gai^x-^ai)) 

'IV- r (8) 

~ g on,x 1(q:P 

^ g-O.SOSVFe 


Hence, the possibility Pq that a man couldn’t find a woman to build up a 
BP can be estimated: 







Figure 4: The fitting of against ^/ns. 


Po= [' f{e)Po{e)de 
Jo 

(9) 

> f f {e){l — 0.808^/ne) de 

Jo 

= 1 - 0.808^ = 0.35, 

which happens to be same with Zhang et al j^. Furthermore, if we substi¬ 
tute Eq. (j5]) into Eq. (|9]), 

g-0.19(VF£)^-0.68(^/F£)+0.02g-0.808^/F 
g-0.19e2-l.488e-l-0.02 ( 10 ) 

= 0 . 60 . 

The result from numerical simulation is Pq = 0.758. 

It is shown in the numerical simulation that in a system consists of n men 
and n women, in which there are 0.758n men who has no potential mate to 
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Figure 5: (a)The amount of BPs Nspijed squares) and (b)the number of men who 
form BPs N^p (blue squares) versus the system scale n in the simulation result. The 
corresponding dashed lines are the linear fits of the data, slop shown in the figure. 

build up a BP, and 0.242u men who have one or more women forming BPs 
with them. In total, there are 0.325n BPs in the system(Fig.5). 

we can think over the probability that more than one woman who form 
BPs with man i: rrii has a possibility of 75.8% to find no one to form a BP 
with him, possibility of 24.2% to have one or more women to build up BPs. 
In the 24.2% of men who is one side of the BPs, if and Wa forms a BP, 
the other n — 1 women besides Wa will have a possibility of 75.8% that could 
not form a. BP with man. That is to say, the possibility that there is one 
and only one woman who builds up a BP with man is 0.242 x 0.758; the 
possibility that there are two or more women who form BPs with man is 
0.242 X 0.242. Among these women, there would be 0.758 x 0.242^ n women 
who forms two BPs. And so on we can see the possibility that a man could 
find exact k women to form k BPs is that: 


Pk = 0.758 X (1 - 0.758)^ A; = 0,1, 2,.... (11) 

Our simulation result is shown in Fig.6. To further ensure with Eq. dm , we 
calculate the total number of BPs: 

CX) 

Nb.p. = 0.758 X 0.242^= xkxn = 0.319n, (12) 

k=l 

which is in good consistency with former numerical simulation. To each man, 
averagely, there would be 0.325 woman who can form a BP with him. 
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Figure 6: The probability of k the number of BPs one has in the systems consist of 
100(blue triangles), 300(purple circles), and 1000(green inverted triangles) men or women. 
The red line is the theoretical prediction by previous analysis. 


In order to get a clearer view of the fore mentioned stability issue, here we 
introduce a parameter which is crucial to the stability of system. Imagine a 
man i who only acquaint only some of the women, for example, a tiny fraction 
k{0 < k 1). In this case, this man i is satished only if he couldn’t hnd a 
BP partner. The possibility which he is satished can be roughly estimated 
as: 


n 

P, = (!-- -^ 1 - 0.325A: ^ 1. (13) 

n 

So the possibility of the whole system being stable is that: 

-I 0.325nfc 

^ 0.72^, (14) 

e 

where K = nk is the average number of women this man acquainted. 

Here we can se^ this parameter k is somehow the key to the stability. As 
mentioned before p , if everybody obeyed the allocation of the matchmaker 
who controls the whole system, this ground state would be the perfect choice 
only considering about the global circumstances. But in fact, there are some 
people of them who do not fully focus on the bigger picture (which is very 
common), then it would be quite possible that someone would like to choose 
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their blocking pair partner rather than the assigned one, should be. If so, 
Considering in the situation that there was no extra constrains like religion, 
moral standards and cultural influences, the ground state matching would be 
broken with several consecutive actions of the pursuit of personal happiness. 
Generally speaking, as we can see from Eq. flT^ . if the marriage society is 
better connected, which means the people will have more relationships with 
heterosexual agents, we can tell that it’s even harder for the system to get 
stable, jsl 

Is there any way leading us to a stable state of the matching problem? 
Yes. Roth et al |l^ proved the existence of such a path. They proved that 
even from arbitrary state, if we keep diminishing the BPs, we will never 
be trapped in a loop, and so we can always hnd a path to the stable state 
anyway. But they didn’t give out the practical method to find such a path, 
it’s very interesting to study how we can get stable state from the ground 
state with this method. 

6. Conclusion 

In this paper we studied numerical simulation result of the ground state 
of Gale-Shapley model with K-M algorithm. With a relaxation on the restric¬ 
tions of initial assigned energy distribution, the simulation result conhrmed 
the theoretical analysis of the ground state using replica method. The un¬ 
solved question about how many blocking pairs are there in the system was 
also shown in the simulations. Furthermore, we discussed the factors which 
affects the stability of the system, and found out that the connectivity of the 
system is crucial to the stability of the ground state. 

In the future, the detailed route path from the erratic ground state, to a 
certain stable state would be of great interests for us to finish. What is the 
micro dynamics of the system when these hopping agents tried to improve 
their own situation and disregarded the cost this bring to the whole system? 
After that, it is also thought-provoking to consider the reality meaning of 
the ground state and the optimal stable solution. May these understanding 
could better ourselves. 
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